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Summary. — We show that there exists a natural Tulczyjew triple in the dynamics 
of objects for which the standard (kinematic) configuration space TM is replaced 
with A"TM. In this framework, which is completely covariant, we derive geometri¬ 
cally phase equations, as well as Euler-Lagrange equations, including nonholonomic 
constraints into the picture. Dynamics of strings and a constrained Plateau problem 
in statics are particular cases of this framework. 

PACS 02.40 - Geometry, differential geometry, topology. 
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1. — Introduction 

This work is a continuation of a research undertaken jointly with W. M. Tulczyjew 
on the Legendre transformation in the dynamics of strings. We show that there exists a 
natural Tulczyjew triple in the dynamics of objects for which the standard (kinematic) 
configuration space TM is replaced with A"TM. To do that, we make use of graded 
bundles of degree n, i.e. objects generalizing vector bundles (for which n = 1). For 
instance, the role of TT*M is played in our approach by the manifold A"TM A" T*M, 
which is canonically a graded bundle of degree n over A"TM. To obtain the dynamics, 
we use the canonical multisymplectic n -\- 1-form on A”T*M, which gives, by the 
contraction, canonical morphism 

Ph'. A”TA"T*M^T* A™T*M. 

The dynamics of strings and the Plateau problem in statics are particular cases of this 
framework. We refer to our work [18] if details of this concept and references are con¬ 
cerned. Here, we add nonholonomic constraints into the picture together with a geometric 
description of the corresponding d’Alembert principle and constrained Euler-Lagrange 
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equation. Sine we are working with Tulczyjew triples, there is no need to introduce 
Poincare-Cartan forms and the presentation of nonholonomic constraints seems to be 
simpler than that known from the literature (c.f. [1, 37], although the idea of d’Alembert 
principle is preserved. The presented approach using multivectors can be viewed also as 
a simplified version of a variational calculus for non-parameterized submanifolds, whose 
full version is technically much more complicated (see e.g. [28]). Some Tulczyjew triples 
for general classical field theories appeared recently in [13, 14, 16]. 

We want to stress that our framework is completely covariant and, what is important, 
not reduced to derivation of just the Euler-Lagrange equations. We present the full pic¬ 
ture, determining the phase space and the phase equations as the meeting point of the 
Lagrange and the Hamilton formalisms, subject to the corresponding Legendre transfor¬ 
mation. The equations are obtained purely geometrically by means of the morphisms 
in the triple from the Lagrangian submanifolds generated by Lagrangians or Hamilto¬ 
nians. In particular, on the Hamiltonian side we do not use, at least explicitly, any 
Poisson brackets. For the basics on Tulczyjew triples and Legendre transformations we 
reccomend [32, 33, 35]. 

Note finally that classical field theory is usually associated with the concept of a 
multisymplectic structure. The multisymplectic approach appeared first in the papers of 
the ‘Polish school’ [8, 23, 24, 34]. Then, it was developed by Gotay, Isennberg, Marsden, 
and others in [11, 12]. The original idea of the multisymplectic structure has been 
thoroughly investigated and developed by many contemporary authors, see e.g. [2, 4, 5, 
6, 7, 31]. The Tulczyjew triple in the context of multisymplectic field theories appeared 
recently in [3] and [27] (see also [38]). A similar picture, however with differences on the 
Hamiltonian side, one can hnd in [9] (see also [10, 26]) and many others; it is not possible 
to list all achievements in this area. 

Note, however, that the multisymplectic structure which appears in this paper is a 
canonical structure on A^T*M and the question about the proper abstract generalization 
of this structure is not discussed here. On the other hand, it would be interesting to 
develop a similar theory for more general objects in the spirit in which mechanics on 
(Lie) algebroids generalizes the classical one for TM (see e.g. [15, 17, 30, 40]). 

2. — The standard Tulczyjew triple 

The canonical symplectic form lom on T*M induces an isomorphism 

/3m : TT*M T*T*M. 

Composing it with 77 .jm, where 


TZe ■■ T*E* T*E 

is the well-known canonical isomorphism o double vector bundles (see e.g. [25, 29, 39]), 
we get the map 


aM ■■ TT*M -)■ T*TM. 

Using the standard coordinates and on TM and T*M, respectively, and 

the adapted coordinates on T*TM and TT*M, we can write 


(1) 


a{x,p,x,p) = {x,x,p,p). 
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This gives rise to the commutative diagram of double vector bundle (iso)morphisms (Tul- 
czyjew triple) 


(2) - TT*M—T*TM 



Note that the mapping um can be obtained directly as the dual to the ‘canonical flip’ 
Km ■ TTM —>■ TTM, which is an isomorphism of two vector bundle structures on TTM: 


(3) 


TTM 


Km 


TTM 



Indeed, the duals of these two vector bundle structures on TTM are T*TM and TT*M, 
and Om can be understood as the dual map of km- 

The map km, as well as um and /3 m, encodes the Lie algebroid structure of TM and 
note that no brackets are needed (cf. [21, 22]). 

3. — The Tulczyjew triple - Lagrangian and Hamiltonian formalisms 

Let now, for a mechanical system, M be the manifold of positions, so that the tangent 
bundle TM represents (kinematic) configurations and T*M is the phase space, and let 
L : TM ^ R be a Lagrangian function. Putting all this into the Tulczyjew triple, we get 
the diagram 



M 


M 


Here, 


VL-.IM rL(x,x) = (x,—) 

ox 
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is the Legendre map, and the submanifold 
(4) V = al,\dL{TM))) 

of TT*M represents the (implicit in general) Lagrange (phase) equation. Recall that, by 
definition, a curve /3 : R —> iV is a solution of an implicit differential equation (differential 
relation) I? C TiV if its tangent prolongation t/3 : M. ^ TN takes values in V. In local 
coordinates. 


^ f/ . .N dL . dL 

p=—, p= — 

so that the phase equations in an implicit form read 




dx 


The Hamiltonian formalism looks analogously. If i? : T*M R is a Hamiltonian 
function, from the Hamiltonian side of the triple 



" M .■: M 


we derive the phase dynamics in the form 

This dynamics is automatically explicit, i.e. generated by the Hamiltonian vector field, 
so it corresponds to a phase dynamics induced by a Lagrangian function only in regular 
cases. In general, one has to use more sophisticated tools like Morse families etc., see 
[36]. In local coordinates, 

. .N . dH . dH 
= <{x,p,x,p)p = - — , ^ = 

so we obtain the standard Hamilton equations. 



4. Nonholonomic constraints and Euler-Lagrange equations 

Let now, 7 : R M be a curve in M (of course, R can be replaced by an open 
interval), and ty : R —>■ TM be its tangent prolongation. It is easy to see that both 
curves, dL o ty and om o t{VL o ty) are curves in T*TM covering ty. Therefore, their 
difference makes sense and, as easily seen, takes values in the annihilator V^TM of 
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the vertical subbundle VTM C TTM. Since V^TM TM Xm T*M, we obtain a 
map 5L^ : M —>■ T*M. This map is usually interpreted as the external force along the 
trajectory. Its value at t € K depends only on the second jet t^ 7 (t) of 7 , so defines the 
variation of the Lagrangian understood as a map 

(5) 5L : T^M T*M , 

where T^M, the second tangent bundle, is the bundle of all second jets of curves 'R ^ M 
at 0 € R. The equation 

(6) SLry = 6L Ot^"f = 0 

is known as the Euler-Lagrange eguation and tells that the curve dT o t 7 corresponds via 
um to an admissible curve in TT*M, i.e. the tangent prolongation of a curve in T*M. 
Here, of course, t ^7 is the second tangent prolongation of 7 to T^M. 

If now, H C TM is an affine subbundle of TM, v(H) is the linear part of A, and 
v(H)° C T*M its annihilator, then we can replace ( 6 ) with 

(7) SL^ G v(A)° , 

which is the d’Alembert principle. It tells that the forces belong to v(yl)°, so make 
no work along the trajectory. The constrained nonholonomic Euler-Lagrange eguations 
associated with the affine nonholonomic constraint represented by A take the form 

(8) tyGH; 

(9) 5L-, G v(H)° . 


In a more traditional form, 
( 10 ) 

( 11 ) 


they can be expressed in local coordinates as 
-a'^))T]l{x{t)) = 0 ; 

dx^ At \dx'^) ’ 


where 77 * = 77 ^( 0 ;) dx"^, i = 1 ,..., dim(v(yl)°), are one-forms generating v(H)°, a{x) = 
{x^,a'^{x)) is an arbitrary section of the affine bundle A —)• M, and Ai(a;) are arbitrary 
coefficients, traditionally interpreted as ’’Lagrange multipliers”. The section a can be 
chosen 0 if the constraints are actually linear. 


5. — The Tulczyjew triple with the kinematic configuration space A"TM 

We want to build a similar framework replacing points with higher dimensional ob¬ 
jects, being motivated by the study of dynamics of one-dimensional non-parameterized 
objects (strings). 

The motion of a system will be given by an n-dimensional submanifold in the manifold 
M (“space-time”). An infinitesimal piece of the motion is the first jet of the submani¬ 
fold. However, this model leads to essential complications even in one-dimensional case 
(relativistic particle). For instance, the infinitesimal action (Lagrangian) is not a func¬ 
tion on first jets, but a section of certain line bundle over the first-jet manifold, a ‘dual’ 
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of the bundle of “first jets with volumes”. Therefore we will take the compromise and 
use for the space of infinitesimal pieces of motions the space of simple n-vectors, which 
represent first jets of n-dimensional submanifolds together with an infinitesimal volume. 
It is technically convenient to extend this space to all n-vectors, i.e. to the vector bundle 
A"TM of n-vectors on M. In this way we get the following principles: 

• A Lagrangian L is a function on infinitesimal motions, L : A"TM —>■ M. If L is pos¬ 
itive homogeneous, the action functional does not depend on the parametrization 
of the submanifold and the corresponding Hamiltonian (if it exists) is a function 
on the dual vector bundle A”T*M (the phase space). 

• The dynamics should be an equation (in general, implicit) for n-dimensional sub¬ 
manifolds in the phase space, i.e. 

V C A^TA”" T*M. 

• A submanifold S in the phase space A"T*M is a solution of V if and only if its 
tangent space TaS” at a G A”T*M is represented by a n-vector from Va- If we use 
a parametrization, then the tangent n-vectors associated with this parametrization 
must belong to V. 

For simplicity, in what follows we will consider the ‘string case’ n = 2, but the 
constructions remain valid for arbitrary n. We will use canonical coordinates 
and (x^jP^i,) on A^TM and A^T*M (with the convention x^^ = — p^i, = —p^fi), 
respectively, representing the decomposition of bivectors: 

A da,. G A^TM, dx>^ A dx'' G A^T^M . 

Since A"T A" T*M is NOT a double vector bundle for n > 1, we start with introducing 
objects naturally generalizing vector bundles (see [19, 20]). 

Definition 1. A graded bundle of degree r is a fibration £ ^ M such that the typical 
fiber is with coordinates . ,y^) which have associated weights (or degrees), 

wi, ... ,Wk = 1 ,...,r, respected by the fiber-bundle change of coordinates. 

We extend weights in fibres, associating weights 0 with basic functions, thus having 
local coordinates {x^,y‘^) consisting of homogeneous functions. Such a structure can be 
conveniently encoded by the weight vector field 

Xs = Way°-dy. , 
a 

whose flow extends to an action R 9 f i-A hj of multiplicative reals 

/^t(x^2/“) = (a:^r“p“). 

In this sense, vector bundles are graded bundles of degree 1, with the Euler vector field 
as the weight vector field. 

Example 2. As canonical examples of graded bundles can serve the higher tangent bun¬ 
dles T^M with the adapted coordinates (a;, x,x,...) of degrees 0,1, 2 ,..., respectively. 
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Example 3. If r : E — M, with affine coordinates is a vector bundle, then 

A^TE is canonically a graded bundle of degree 2 with respect to the projection 

A^Tt : A^TE -A A^TM . 


The adapted coordinates ,y °‘= —z'^^, coming from 
the decomposition of a bivector 


A^TTM Au= 


d d 

dx^ ^ dx’^ ^ 


ah 


d d 1 cd ^ '9 

dx'^ dy^ 2 dy'^ y'^ 


are of degrees 0,1, 0,1, 2, respectively. 

Definition 4. Like for double vector bundles, two structures of a graded bundle are 
called compatible, i.e. define a double graded bundle, if the corresponding weight vector 
fields commute. 

Example 5. Our canonical example will be A^TE with the diagram of bundles and their 
morphisms 



In particular, A^T A^ T*M is a double graded-vector bundle 


A^T A^ J*M 




A^JM 



On A^T*M, we have the canonical Liouville 2-form: 

^li = A dx‘', 

inducing the canonical multisymplectic form 

^ -d dx>^ A da;*". 

The multisymplectic form, via the contraction, induces in turn the double graded bundle 
morphism 

Ph : a2 T a2 T*M ^ T* a2 T*M 
: u i-A i„a;|f. 
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In local coordinates (summation convention used), 

,yep,P'rSec) = ,P\k,■ 

Using now the canonical isomorphism of double vector bundles 
7^ = 7^^2TM : T* T*M ^ T* TM , 
we can define aj^ = TZ which is another double graded bundle morphism, 

alt'. a2Ta2T*M^T*a2TM, 

of double graded bundles over A^TM and A^T*M. In local coordinates, 
o^hix^,P\K,x'''^,yep,P7Sec) = {x^,x'''^,y'^p,p\K.) ■ 

The map aj^ can also be obtained as a certain ‘dual’ of the canonical isomorphism 

: T A^ TM A^TTM , 

generalizing the canonical flip (3). Combining the maps and we get the following 
Tulczyjew triple for multivector bundles, consistsing of double graded bundle morphisms: 

T* a2 T*M A^T A^ T*MT* A^ TM 



We have a straightforward generalization for all integer n > 1 replacing 2: 

Ph - 


T* A" T*M- 


■ A”TA”T*M- 


M 


■ T* A” TM 



The map 


: A" T A” T*M ^ T* A" T*M , 
: u i-A IuUjIp , 
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comes from the canonical multisymplectic (n + l)-form on A"T*M, being 

the differential of the canonical Liouville n-form 

9 m = A • • • A dx'"" . 

The map is just the composition of /Sjjf with the canonical isomorphism of double 
vector bundles T* A” T*M and T* A" TM. 

6. String dynamics 

The way of obtaining the implicit phase dynamics V as a submanifold of A^T A^T*M, 
from a Lagrangian L : A^TM —>■ R or from a Hamiltonian H : A^T*M —>■ R, is now 
fully analogous to the case of the standard Tulczyjew triple. For the Lagrangian, the 
corresponding diagram reads 


V A^T A^ T*M - — -^ T* a 2 TM dL 



and 

In local coordinates, 

P= |(x'^,PAK,x"",y^^,:P^*c) : y^p = PA. = • 

On the Hamiltonian side, we get 



and 


V = {f]l,)-\dH{A^T*M)). 
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In local coordinates, 


r 3M 

V = : Vvp = 


dH 

^Pva 


that explicitly produces the Hamilton equations. Again, this framework can be extended 
to Morse families. 


7. — The nonholonomic Euler-Lagrange equations for strings 

To define nonholonomic Euler-Lagrange equations for strings, consider a surface 

S' : 9 {t, s) 1-9 s)) 

in M and its bi-tangent prolongation 

A^tS : ^ A^TAf, A^tS = HS A t,S, 

where tt and t^ are tangent prolongations of curves with respect to variables t and s, 
respectively. 

It is easy to see that both parameterized surfaces, dLoA^tS and a^oA^t('PLoA^tS) 
in T* A^ TM cover A^tS. Therefore, their difference makes sense and, as easily seen, 
takes values in the annihilator IA° A^ TM of the vertical subbundle V A^TM C T A^ T M. 
Since A^ TM ~ A^TM Xm T*M, we obtain a map SLs : -A T*M. The above 

map is interpreted as external forces along the string trajectory S. Its value at (t, s) 
depends on the second jet s) of S only, so defines the variation of the Lagrangian 

understood as a map 

(12) <5L ; J2(R2^M) ^T*M, 

where Jg(R^,M) is the bundle of all second jets of maps —)• M at 0 G R^. The 

equation 

(13) SLs = 0 

we will call the Euler-Lagrange equation. It tells that the surface dL o A^tS corresponds 
via to an admissible surface in A^T A^ T*M, i.e. the bi-tangent prolongation of a 
parameterized surface in A^T*M. 

If now, A C A^TM is an affine subbundle of A^TM, v(A) is the linear part of A, and 
v(A)° C T*M its annihilator defined as the set 

(14) v(A)° = {r/^ G T*M | ir,^u^ = 0 for all G v(A) 2 ,} , 
then we can replace (6) with 

(15) dLs€v(A)°, 

which is the dAlembert principle. It tells that the forces SLs belong to v(A)°, so make 
no work along the trajectory. Indeed, that a one form rj vanishes on the surface S is 
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equivalent to the fact that its contraction with X A F is zero, for any vector fields X, Y, 
tangent to S. 

In other words, S satisfies the Euler-Lagrange equations if the image by dL of its 
prolongation to A^TM, 


{t,s) i-A 




dx^ dx^ dx^ dx^ \ 

dt ds ds dt ) ' 


is a^-related to the prolongation of the surface VL o A^TS', living in the phase space 
T*M, to A^T A^ T*M. 

The constrained Euler-Lagrange equation associated with the affine nonholonomic 
constraint represented by A take the form 


(16) 

(17) 


A^tS G v(A); 
5Ls G v(A)° . 


In coordinates, they read 


(18) 

(19) 


fdS^^ dS'' 

dS'^ dS^^ 


\ dt 

ds 

dt ds 

dL 

dS^^ 

d f dL \ 

dS^^ d 

dx’' 

dt 

ds V dx^^’^) 

^ ds dt 


S) =0; 
dL 


dxf^’' 


= Kvl , 


where, as before, 17 ® = ry^cc) dx'^, i = 1,..., dim(v(A)°), are one forms generating v(Al)°, 
a = {x°^a^'^ix)) is an arbitrary section of the affine bundle A —>• M, and Xi{x) are 
arbitrary coefficient functions. 

Example 6. In the dynamics of strings, the manifold of infinitesimal configurations is 
A^TM, where M is the space time with the Lorentz metric g. This metric induces a 
scalar product h in fibers of A^TM, so that for 


d d 
dxf^ dx’' 


d d 
dxf^ dx’^ 


we have 


(u|w) = x"^^ , 


where 


h 




The Lagrangian is a function of the volume with respect to this metric, the so called 
Nambu-Goto Lagrangian, 


L{w) = \J {w\w) 




which is defined on the open submanifold of positive bivectors. 
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The dynamics V C A^T T*M is the inverse image by a\j of the lagrangian sub¬ 
manifold dL(A^TM) and it is described by the Lagrange (phase) equations 


where 


Vau 


1 ■fJ.K 'Act 

2p dx‘' > 




Xk, 


The dynamics V is also the inverse image by of the lagrangian submanifold in T* A^ 
T*M, generated by the Morse family 

H : a2T*MxM+ 

: (p,r) i-A r{^(p\p) - 1 ) . 

Example 7. In the case of a minimal surface, i.e. the Plateau problem^ we replace the 
Lorentz metric with a positively defined one. 

In particular, if M = = {(x^ = x,x^ = = z)} with the Euclidean metric, the 

Lagrangian reads 


L{x>^,x^^) 



The Euler-Lagrange equation for surfaces 


( 20 ) 


S : {x,y) i-A {x,y,z{x,y)) 


provides the well-known equation for minimal surfaces, found already by Lagrange, 


( 21 ) 



A 

dy 



= 0 . 


In another form, 

(1 - 1 - Z^')Zyy ‘IZxZyZxy -^ (1 Zy^ZxX — 0 • 

Example 8. For the Plateau problem as above, let us put nonholonomic constraints 
given by the affine bundle 


A = {^x^^y + f{x,y,z){dx - dy) Adz\f £ (17°“ (K^)} . 

It can be easily seen that v(^)° is generated by the one-form 77 = dx-|-d?/. It is the matter 
of direct computations to show that the surface (20) satisfies (18) and (19) if and only if 
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(21) is satisfied together with Zx, = Zy. But the latter means that z{x, y) = F{x + y) for 
some function F of one variable, and we end up with the equation 



which is equivalent to F" = 0. Hence, the solutions of our nonholonomic string equation 
are just planes being the graphs of linear functions u{x, y) = a{x + y) + b. 
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